format 8 

* Functions & variables part:
var laF la/xFO
var laC (1-xFO)*laF/xCF
var laP laC*xCD
var laD laC*xCD
var WF 1 - (4*(laF/muF)^4/(24*(4-(laF/muF))))*(1+(laF/muF)+0.5*((laF/muF)^2)+(1/6)*((laF/muF)^3)+(1/24)*((laF/muF)^4)*((4*muF)/(4*muF-laF)))^-1
var WP 1-((laP/muP)^2/(1-(laP/muP)))*(1+(laP/muP)+0.5*((laP/muP)^2)*((2*muP)/(2*muP-laP)))^-1

bind 
muT   0.2
muP   0.05
muF   0.375
muC   1
muD   2
CF   4
CP   2
xFO   0.5
xFC   0.5
xCF   0.46
xCD   0.33
xCP   0.21
la   1/30
end

markov responseblue 
T F1 muT
F1 Done WF*xFO*muF
F1 F2 muF*(1-WF)
F1 C xFC*muF*WF
C P1 xCP*(muC-laC)
C F1 xCF*(muC-laC)
C D xCD*(muC-laC)
P1 C WP*muP
P1 P2 muP*(1-WP)
D C (muD-laD)
F2 Done xFO*(CF*muF-laF)
F2 C xFC*(CF*muF-laF)
P2 C (CP*muP-laP)
end
T 1
end

expr laF
expr WF
expr laP
expr WP
expr laC
expr laD

end

